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To quantify the effect of decoherence in quantum measure- 
ments, it is desirable to measure not merely the square modu- 
lus of the spatial wavefunction, but the entire density matrix, 
whose phases carry information about momentum and how 
pure the state is. An experimental setup is presented which 
can measure the density matrix (or equivalently, the Wigner 
function) of a beam of identically prepared charged particles 
to an arbitrary accuracy, limited only by count statistics and 
detector resolution. The particles enter into an electric field 
causing simple harmonic oscillation in the transverse direc- 
tion. This corresponds to rotating the Wigner function in 
phase space. With a slidable detector, the marginal distribu- 
tion of the Wigner function can be measured from all angles. 
Thus the phase-space tomography formalism can be used to 
recover the Wigner function by the standard inversion of the 
Radon transform. By applying this technique to for instance 
double-slit experiments with various degrees of environment- 
induced decoherence, it should be possible to make our under- 
standing of decoherence and apparent wave-function collapse 
less qualitative and more quantitative. 



I. INTRODUCTION 

The problem of how to interpret measurement in quan- 
tum mechanics has caused intense debate ever since 1925, 
and shows little sign of abating. However, driven by 
experimental progress in for instance low-temperature 
physics and quantum optics, the debate is changing in 
character, becoming more quantitative than qualitative. 
The perennial question of whether the wavefunction of 
some given system evolves according to the Schrodinger 
equation or for all practical purposes collapses need no 
longer be discarded as mere metaphysics. Rather, it can 
often be answered experimentally (e.g. jl]|2)), and in some 
cases even answered by direct computations of the impact 
of the environment upon the system (e.g. |^,||]), quanti- 
fying the apparent wave-function collapse known as de- 
coherence f5p|^]. 

Our knowledge of the state of a quantum system is 
completely described by its density matrix p |9| . It gener- 
alizes the wavefunction description by incorporating our 
lack of knowledge as to what pure state the system is 
actually in. To be be able to further refine our under- 
standing of the measurement process, decoherence, etc., 
it is clearly desirable to be able to accurately measure 
this key quantity p, and several formal methods have 
been proposed for doing this fl0| - |l4) . Our apparatus is 



based on the technique known as "phase-space tomogra- 
phy" |l5|-|l8| (also rediscoved independently by the au- 
thor) which has been successfully applied to a number of 
cases involving measurements of the electromagnetic field 
|l9|-p4| . The purpose of this paper is to show how phase- 
space tomography can be applied to one of the most basic 
cases in quantum mechanics: the spatial density matrix 
of a charged particle. 

II. THE APPARATUS 

The apparatus is illustrated in Figure [l]. It consists 
of a slidable particle detector inside of a shielded box 
where an electric field makes the entering charged par- 
ticles (which we will take to be electrons, for definite- 
ness) feel a simple harmonic oscillator potential in the 
x-direction. Inside the box, the Coulomb potential is 

<f> = ^ 2 (y 2 -- 2 )- W 

The box consists of a large number of metal plates, insu- 
lated from one another, and since V 2 (f> — 0, this desired 
field configuration is readily arranged by fixing the po- 
tentials of these plates at the appropriate values as shown 
in the figure. 



The Wigner phase space distribution W of a ID quan- 
tum particle is related to its density matrix p by |25[| 

W(x,p) = ±Jp(x-~,x+±)e* u du, (2) 

i.e., essentially by an inverse Fourier transform in the off- 
diagonal direction followed by a 45° rotation. Since the 
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FIG. 1. The density matrix measurement apparatus. 
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Hamiltonian for an electron inside our box is quadratic 
in positions and momenta, it is well-known [§5-27 that 
the time-evolution of the phase space distribution W is 
purely classical, given by the Liouville equation. 

The motion is clearly independent in the x—, y— and 
z— directions, corresponding to a harmonic oscillator, an 
upside down harmonic oscillator and a free particle, re- 
spectively. Since slower electrons will curve more in the 
magnetic field to the left of the box, it is easy to arrange 
for our electron beam to be highly monochromatic, with 
Ap y <C ip y )- In this limit, the marginal Wigner distri- 
bution for the a;— direction will evolve as 

W t (x,p) — Wo(xcos9 — p sin 8, x sin 9 +pcos9), (3) 

where we have defined 9 = cut and 



'2V \q e 



L 2 7 



(4) 



We have chosen our units so that mu> = 1, to avoids 
cumbersome conversion factors between position and mo- 
mentum. In other words, the time evolution simply cor- 
responds to a clockwise rotation of the Wigner function, 
as shown in Figure ^. 



respectively. In X-ray tomography, one measures the in- 
tegral of the integrated X-ray opacity through say a pa- 
tient's head, as seen from a large number of angles 9, 
and then wishes to reconstruct the 2D cross section. We 
can do "phase space tomography" and obtain "X-ray im- 
ages" of the Wigner function from different angles by 
sliding the detector (changing 9). 





III. HOW IT WORKS 

Defining t — as the time when a particle passes y = 
(the origin of coordinates is at the center of the box) and 
p y o as (p y ) at that time, we can make the identification 



y = p y0 sinh( 



(5) 



since Ap y <C (p y ) at all times (particles with small y- 
momentum have been rejected by the collimators). When 
the slidable detector is positioned at y, it will thus detect 
the probability density 



p(x) 



W t (x,p)dp, 



(6) 



where t = 9/cu is given by equation (|5|). Substituting 
equation ([}]) into equation (||), we can rewrite the integral 

as 



p[x) = 



<5(n-r-x)W (r), 



(7) 



where r = (x,p) and n = (cos#, sin#). We recognize 
the right hand side of equation (]?]) as the very definition 
of the Radon transform of Wo, conventionally denoted 
Wo(ri,x). The Radon transform has a simple geometri- 
cal interpretation: Wo(n,x) is the marginal distribution 
of Wo projected onto a line parallel to the vector n. The 
two "shadows" shown in Figure || are the marginal distri- 
butions for x and p, corresponding to 9 = and 9 = tt/2, 



FIG. 2. How "phase-space tomography" works. Shifting 
the detector to the right corresponds to rotating the Wigner 
function as the arrow shows, so that the detector measures 
the marginal distribution seen from a different angle in phase 
space. The marginal position and momentum distributions 
(shaded) correspond to 9 — and 8 — n/2, respectively. 

The Radon transform is closely related to the Fourier 
transform, and it can we shown that |2q| 



e lrx W (ri,x)dx = W (nr), 



where Wq is the 2D Fourier transform of Wo ■ The Fourier 
transform of the Wigner function is often called the char- 
acteristic function, and substituting equation (0) into 
equation (||), we thus find that the characteristic func- 
tion is given by simply 



Wo(Kr) 



c p(x)dx=p*{r). 



(9) 



In other words, when we Fourier transform the proba- 
bility distribution measured by the detector, we obtain 
a radial strip of the characteristic function. Sliding the 
detector to a new location gives another radial strip, etc. 
When we have covered phase space densely enough with 
such strips, we perform a 2D inverse Fourier transform 
(with respect to x and p this time, not with respect to 
the radius r = \J x 2 + p 2 ) and obtain our desired Wigner 
function. Alternatively, we can compute the density ma- 
trix directly from W as 
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p(x, x') = 7-— / W(k, x' - x)e l(x+x '^ k / 2 dk. (10) 
2tt J 

Figure ^ shows the Wigner function corresponding to 
a pure state where the wavefunction is the sum of two 
Gaussians separated by 10 standard deviations, a very 
crude model of the state of an electron after passing 
through a double slit. When the detector is at the center 
of the box, at y = & = 0, it would measure the double- 
humped marginal distribution for x shown. Moving to 
the right in the box, the Wigner function rotates, the two 
humps start overlapping, and interference fringes begin 
to appear on the detector, finally giving the marginal dis- 
tribution for p (also shown in Figure ^ when 9 = tt/2. 
Destroying the coherence of the electrons (the purity of 
the quantum state) would correspond to removing the 
oscillatory center of the Wigner function in the figure || . 
Thus no fringes would appear as we moved the detector, 
and the two Gaussians would merely add incoherently 
when they overlapped at 9 = tt/2. 

IV. REAL WORLD ISSUES 
A. How to chose the voltage 



This means that its Fourier transform, the Wigner func- 
tion, will have a resolution of order Aa;. (Recall that the 
x — p conversion factor is m e oj, so the momentum res- 
olution Ap = m e ujAx.) Features on smaller scales will 
be washed out, and it is well-known that smearing the 
Wigner function tends to decrease the purity of a state, 
much in the same way as decoherence does. In other 
words, to be able to measure quantum effects with our 
device, not merely classical-looking mixed states, it is 
crucial that interference fringes be present on a scale ex- 
ceeding our resolution. We can arrange this in a va- 
riety of ways, corresponding to various known ways of 
demonstrating visible electron interference patterns, and 
placing the corresponding contraptions to the left of the 
harmonic oscillator box. We might for instance place 
a crystal near the entrance of the box, whose electron 
diffraction pattern could be readily detectable. Alterna- 
tively, we could replace it by a microfabricated double 
slit magnified by electrostatic cylinder lenses, as is done 
in the famous electron version of Young's double-slit ex- 
periment p9[ . Also, the single opening in the box can of 
course be replaced by many well-separated openings, as 
long as they are small enough to leave the interior field 
approximately of the form given in equation (pi) . 



Let us rewrite equation (|J) as 
/ sinh 9 



V = 



\ sinh 9„ 



L, 



where 9 r . 



sinh- 1 [(2V |g e 



m , 



a/2 



/P. 



y0\ 



(11) 



To be able to 



"X-ray" the Wigner function from all angles — tt/2 < 9 < 
tt/2, we clearly want 9 max > tt/2. On the other hand, it is 
not feasible to make 9 ma x much larger than this, since the 
required voltage Vq eventually grows exponentially with 
9 max- We therefore suggest chosing 9 max only slightly 
larger than tt/2, to allow for the finite thickness off the 
detector and room for optional double slits etc. on the 
left hand side. This produces classical trajectories such 



c°Scos sinh 



2. 3 |; 



(12) 



the beam curve in Figure [l], and means that Vq will be 
about five times the potential that initially accelerated 
the incident electrons. 



B. Overcoming the resolution 



C. Other constraints 

The main additional constraint is that we must ensure 
that the Schrodinger time evolution of equation (|j|) is in- 
deed valid inside the box. Firstly, this clearly requires 
that the electrons evolve as an isolated system, for in- 
stance that the vacuum in the box be hard enough that 
the effect of air molecules can be neglected. Secondly, the 
electrons must only "feel" the potential of equation (Q), 
i.e., stay inside the box. A problem of this type would of 
course immediately be noticed, as counts near the edge 
of the detector. 

What number n of different detector locations should 
we use? We saw that W is accurately measured out to a 
radius 1 /Ax. Since each radial strip is a Fourier trans- 
form p, the radial resolution will be limited by the inverse 
length of the detector, 1/2L. Since the resolution in the 
angular direction at the radius 1/Ax is just 2tt /nAx, and 
for economy we want the two resolutions to be roughly 
equal, we should choose n and Aa; such that n ~ L/Ax, 
i.e., so that the number of y- values at which we measure 
roughly equals the number of resolution elements on the 
detector. 



If the detector registers N hits, we must smooth the 
measured p on the scale of their mean separation 2L/N 
to suppress Poisson shot noise. We thus define the spa- 
tial resolution Ace as either 2L/N or the intrinsic reso- 
lution of the detector, whichever is larger. The W that 
we measure will thus be near the true W within a cir- 
cle of radius ~ \/Ax but tend to zero at larger radii. 



V. DISCUSSION 

We have presented a method for measuring the ID spa- 
tial density matrix (equivalently, the Wigner function) of 
a beam of identically prepared charged particles. Specifi- 
cally, we measure the Wigner function that describes the 
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ensemble of particles when they are at y = 0, half way 
through the box. Some clarifications are in order here: 

• The Wigner function at another y-coordinate is 
obtained by simply rotating the measured Wigner 
function by the appropriate amount. 

• The reader may feel uneasy about the fact that our 
measurement technique assumes the validity of the 
Schrodinger equation. However, this per se is not 
that different from say measuring the velocity of a 
classical object, which requires position measure- 
ments at two different times and the assumption 
that Newton's law of motion is valid during the in- 
terval. 

• By the very definition of the density matrix, we can 
never measure it for a single particle, merely for an 
ensemble. 

• The condition "identically prepared" is in a sense 
fulfilled by definition: if the particles in the ensem- 
ble are in fact not all in the same state, but we are 
unaware of this, this lack of knowledge will merely 
be reflected in the density matrix we measure. 

The implementation of phase-space tomography pre- 
sented here can obviously be generalized in a number 
of ways. For instance, it may be feasible to measure p 
for neutral particles (discussed by p3pl| ]) by replacing 
the box by a Stern-Gerlach type apparatus coupling to 
the spin, or by an electric field whose gradient couples to 
the dipole moment of the particles, in such a way as to 
produce harmonic motion in the x-direction. Indeed, it is 
easy to show that our Radon transform approach can be 
applied even without any external force field, for free par- 
ticle time evolution. In this case, however, we never ob- 
tain quite all radial strips of W , since the Wigner function 
will not rotate but shear, and 9 = ±ir/2 will correspond 
to t = ±00. More general stationary and time-varying 
potentials can of course also be used — the advantage 
of our quadratic potential was merely that the resulting 
inversion problem was linear and easy to solve. 

As to the types of beams for which p can be measured, 
the variations are of course many as well, employing 
various combinations of the above-mentioned diffracting 
crystals, double slits, electrostatic and magnetic lenses, 
etc., and adding various sources of decoherence. 

In summary, it is hoped that this rather versatile tech- 
nique can help us continue the trend mentioned in the 
introduction, and make our understanding of decoher- 
ence, quantum measurement and apparent wavefunction 
collapse less qualitative and more quantitative. 
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N.B. The latest version of this paper and related work 
is available from 

h t t p://www. sns.ias.edu/' max /radon. html 
(faster from the US) and from 

h t t p://www.mpa- garching.mpg.de/' max/radon. html 
(faster from Europe). Note that figure 1 will print in 
color if your printer supports it. 
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